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Abstract 

The c- vortex ensembles are constructed by means of the recently proposed 
cooling method which gradually removes the SU{2)/Z2 coset fields from the 
SU{2) lattice configurations and which thus reveals the Z2 vortex vacuum 
texture. Using Teper's blocking method, the screening masses of the 0^ and 
the 2"'" glueball is calculated from these vortex ensembles and compared with 
the masses obtained from full configurations. The masses of either case agree 
within the achieved numerical accuracy of 10%. As a byproduct, we find 
that the overlaps of the Teper operators with the glueball wavefunctions are 
significantly larger in the case of the c- vortex ensembles. 
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Introduction. The idea that the center part of the SU{N) gauge configurations 
are important for the confinement of quarks dates back to the late seventies []T[. 
Indeed, the proposal that the vortex free energy serves as an order parameter for 
quark (de-) confinement was recently confirmed in the case of a SU{2) gauge group 
by a large scale lattice simulation g]. Subsequently, Mach and collaborators explic- 
itly constructed a vortex signature of the Yang-Mills vacuum from gauge invariant 
variables 0|. It was observed at that time that random fiuctuations of the vortex 
structure disorders the Wilson loop and, hence, provides quark confinement. With 
the increase of the computer performance in the mid eighties, many research ef- 
forts were devoted to substantiate this idea on a quantitative level 0, || . One finds 
that projecting the full lattice configurations onto its vortex content reproduces the 
string tension of the static quark anti-quark potential. Recently, it was pointed out 
that not only the asymptotic behavior of the potential but also the short range part 
which is due to gluon exchange is reproduced if the plaquette is used for a definition 
of the vortex ensemble [Q. Moreover, the properties of the vortices arising from 
the plaquette projection technique strongly depend on the size of the lattice spacing 
thus rendering a continuum interpretation of the latter vortices cumbersome [Q. 

A significant upturn of the vortex picture of quark confinement occurred with the 
construction of the p-vortices which are defined after adopting the so-called center 
gauge [0] by projecting the gauge fixed link variables onto center elements p|, 0. 
The fact that Yang- Mills lattice configurations which were reduced to their p-vortex 
content reproduce the string tension is often referred to as center dominance of the 
string tension. Moreover, it was observed that the p-vortices are sensible degrees 
of freedom in the continuum limit p, |^: the (area) density of the p- vortices as 
well as their binary interactions extrapolate to the continuum. The p-vortex picture 
of the Yang-Mills ground state also provides an appealing explanation of the de- 
confinement phase transition at finite temperatures |jTO|. 



It was pointed out that the center gauge fixing which is prior to identify the physical 
vortex structure might be plagued by a so-called practical Gribov problem |TT 



In addition, it was observed that the vortex properties are quite sensitive to the 
finite size of the lattice volume [Q. For avoiding the practical Gribov and related 
problems, a gauge invariant definition of the vortex vacuum texture was achieved 
by employing a new self-restricted cooling procedure which diminishes the coset 



fields while leaving the center degrees of freedom un-changed |13]. For rating the 
phenomenological importance of these, say, c- vortices, center dominance of the string 
tension was verified. In addition, the SU{2) action density which is carried by the 
c-vortex vacuum texture properly extrapolates to the continuum limit and, hence, 
gives rise to a mass dimension four condensate which features in the operator product 
expansion 
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In this letter, we will investigate the question whether the c-vortex vacuum texture, 
besides the appealing picture for quark confinement and their hypothetical signature 



in high energy hadron collision experiments ||13[ , also provides a quantitative picture 
of the low lying excitations of pure Yang-Mills theory. For these purposes, we 
will calculate for the first time the correlation function for the O"*" and 2"*" glueball 
channel, respectively, using vortex projected configurations. For this investigation, 
we confine ourselves to the most simple, but academic case of a pure SU{2) gauge 
theory. An analogous investigation was recently performed using ablian projection 
to the maximum ablian gauge in this case, it was observed that abelian 



or 



even monopole) projected configurations still reproduce the 0"'"/2+glueball masses 
known from the full theory. 

Glueball correlation functions. Glueball screening masses are calculated 
from the correlation functions 



by analyzing the exponential decrease of C{t) at asymptotic values of t, i.e. 



C{t) oc exp|— mgt| for t ^ l/rUg 



(2) 



Thereby <f){x) is a combination of the link variables which carry the quantum num- 
bers of the glueball under investigation. For choosing a function which gener- 
ates sufficient overlap with the glueball state, we closely follow the pioneering work 
of Teper [jl5i and define composite link variables 



r(^-i), 



X + i) 



(3) 



{x)Ul'''^\x + k)Ul''^'\x + k + i)U, 



X + i + i] 



where i 



ensure U- U- 



1 . . .3,k G —3, —2, —1, 1, 2, 3 and where A/" is a normalization factor to 
1. The link variables U^^\x) are defined on a coarser lattice of 
size {N^^'^ /2y^Nt where N^'^^ and A^^^ are the number of lattice points in the spatial 
directions and in the time direction, respectively, of the finer lattice. The level 
A^ = corresponds to the finest level of a A^j^-* = A^,. The glueball operators are 
defined by means of the plaquette 
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tr 
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(4) 
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Figure 1: Renormalized correlator CR{t) (^) in the O"*" and 2^ glueball channel, 
respectively, obtained from full SU (2) configurations. 



Here, we use the "blocking" level N = 2 throughout the letter. Using the com- 
posite link method for the glueball operators, the glueball screening masses were 
successfully obtained even for the more realistic case of the gauge group SU{3) ]TB| . 

In the present paper, we use a A^^ = 8, A^^ = 16 lattice and the standard heat 
bath algorithm to generate the link configurations according to a probability dis- 
tribution provided by the Wilson action. Using (3 = A/ (where g is the bare 
gauge coupling) up to values 2.5, we are aware that finite size effects become vis- 
ible. Our point is, however, to compare the glueball masses calculated from full 
configurations with those which were obtained by reducing the lattice variables to 
vortex ensembles rather than to perform a high precision extrapolation to the in- 
finite volume limit. 15000 measurements separated by 10 Monte-Carlo sweeps to 
reduce auto-correlations were performed to estimate C{t) (0). In order to express 
t in physical units, we use the /3-dependence of the lattice spacing a predicted by 
one-loop perturbation theory, i.e. 

aa\(5) = 0.12exp|-^(/3-2.3)|, (6) 

where the string tension a = (440MeV)^ was used as reference scale. 

Since C(0) is the expectation value of a composite field, C(0) acquires additional 
divergencies even if C(t 7^ 0) is renormalized finite |]T^. Hence, we refrain from 
normalizing C(0) to 1, but demand 

CrU) = ZlCit), Cnito) = 1, (7) 
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where to is the renormahzation point. From our numerical simulation, the function 
C(t) at the points t = a, . . . , 4a obeys the exponential law to good accuracy. 
This fact allows to evaluate C(to) by interpolation where to = 0.2 fm is used through- 
out this paper. The final result Cj?(t) is shown in figure |I| for j3 = 2.3, 2.4, 2.5. The 
corresponding data points are satisfactorily close to a single exponential curve, thus 
establishing a renormahzation group invariant screening mass. Note that Cr{0) 
changes if different /? values are used, thus reflecting the additional divergency as- 
sociated with the composite operator. The straight lines in figure (|lD represent 
exponential fits to C/j(t) for each f3. Averaging over the screening masses obtained 
by the fit for a given (3, we find 

mo+ ^ 1.67 ± 0.11 GeV , m2+ ^ 2.30 ± 0.08 GeV . (8) 

The uncertainties provided in (j^) comprise statistical as well as systematic errors 
due to the extrapolation to the continuum limit. The masses (|D are consistent with 
the data presented in O, IP 



Glueball masses from the c-vortex ensembles. In order to reveal the vortex 
vacuum structure, we employ the self-restricted cooling procedure proposed in |TB 



and fractionize the gauge group SU{2)=Z2 x S0{3). The corresponding degrees 
of freedom are center vortices and coset fields. The coset part of the SU{2) link 
variables U^{x) is isomorphic to the adjoint link 

Ofix) = ^tr|t/,(x) Ulix) r"} = O^^f^] , Of{x) G ^0(3) , (9) 

which can be uniquely represented by a gauge vector potential Afj^{x) in the standard 
fashion. For removing gluonic (coset) degrees of freedom from SU{2) configurations, 
the gluonic action density per link is defined by |T3[ 



= E |i - \tTAO,,ix)] = \y: p^i^] + <^i^') ^ (10) 

where 0^^{x) is the plaquette calculated in terms of the 5*0(3) link elements 0^(x) 
(^. The sum over u runs from —4. . .4. is the (continuum) field strength 

functional of the (continuum) gluon fields A^{x). 

Cooling is performed by locally reducing the total gluonic action ( p!OD with respect to 
the fields 0^{x). Further cooling of the adjoint link 0^{x) is rejected if the gluonic 
action is smaller than some threshold value 

sj(x) < SK^aK (11) 

Thereby k is a gauge invariant cooling scale of mass dimension one. The cooling 
procedure stops if the gluonic action density ([T0|) locally has dropped below the 
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Figure 2: Renormalized correlator C/j(t) (^) in the O"*" and 2"^ glueball channel, 
respectively, obtained from c-vortex ensembles. 



critical value specified by k. Details for the practical application of the cooling 
procedure can be found in ||13[. For k = 0, the cooling procedure completely removes 
the gluon fields from the SU(2) lattice configurations leaving only gauge equivalents 
of 0^{x) = 1. In fact, even for ^ a the SU{2) action density is largely clustered 
along 2-dimensional vortex world sheets, and the short range force between a static 
quark anti-quark pair is strongly affected. This is expected since the behavior at 
small distances is dominated by the exchange of gluons, which are already partially 
eliminated by cooling. Throughout this paper, the configurations which emerge after 
"gluon" cooling with k,"^ < a are labeled c-vortex ensembles. 

In order to investigate the c-vortex dominance of the glueball screening masses, we 
evaluate the glueball field combinations (1) and (J^) using c-vortex ensembles ob- 
tained after adjoint cooling. Note that the action density of the gluon (coset) fields 
is by construction limited to < (440 MeV)'^ (see (0) and (|ll|)) implying that the 
large SU{2) action density carried by the c-vortices is required to sustain screening 
masses of order 1.5 GeV. We employ 5000 measurements to obtained the renormal- 
ized correlation function Cnit) calculated with adjoint cooled configurations, using 
k"^ = a and k"^ = 0.5 a, respectively. The result is shown in figure ^ Again, the 
results for several values of f3 are consistent with a single exponential law (§) refiect- 
ing proper scaling towards the continuum limit. The straight line shown in figure ^ 
is the fit to the data obtained from full configurations for /3 = 2.4 (see figure 
We therefore find that the screening masses for the 0"'"/2"'" glueball calculated from 
full configurations and c-vortex ensembles, respectively, coincide within the achieved 
numerical accuracy. 
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Note that the coohng procedure strongly effects the value Cr{0). This is expected 
since the cooling procedure eliminates UV-divergencies which are generated by gluon 
(coset) fields of high (5*0(3)) action density, and therefore alleviates the divergen- 
cies of the composite operator. Furthermore, we point out that roughly the same 
statistical error of the screening masses was achieved in the case of the c-vortex 
ensembles with a number of measurements which is a factor of three less than the 
number of measurements employed in the case of full SU{2) configurations. This 
improvement is due to an enhanced overlap of the "wavefunctions" (^ and (|]) with 



the glueball wavefunctions once the adjoint cooling operates. Following Teper |15 
for an estimate of the overlap, we find 



C{2a) 



C(0) 



for the 0+ glueball, and 



C(2a) 



4.6 



C{2a) 



C(0) 



13.9 



C(0) 



C(2a) 



P = 2.4, 



K 



a 



full 



(12) 



C(0) 



full 



/3 



2.4, 



a 



(13) 



for the 2+ glueball, respectively. 



Conclusions. Recently, a self-restricted cooling method was proposed |1^ which 
gradually removes the SU(2)/Z2 gluon (coset) fields from the SU{2) lattice con- 
figurations paving the way to gauge invariant (c-)vortex ensembles. Self-restriction 
ensures that the local gluonic action density does not exceed the cooling scale k. 
By definition, c-vortex configurations are obtained for the choice k < \fa, where 
a = (440 MeV)^ is the string tension. In the present letter, we have studied the 



screening masses for the O"*" and 2"*" glueballs extracted from Teper correlators |15 
^/o' and K 



Using K = and n = 0.7 y/a, we have shown that these screening masses are 
insensitive to the new cooling method hence providing evidence that the 0"'' and 2+ 
glueball masses are dominated by c-vortex configurations. 
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